Empirical estimation the effects of flow on thermal photon angular distribution and 
spectra in nucleus-nucleus collisions at RHIC and LHC 
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Experimental data for hadron radial and elliptic flow are used to investigate their influence on 
the shape of the thermal photon spectrum at RHIC. Leaving alone the actual mechanism of photon 
production and its time evolution, we concentrate on the spectrum shape. Radial and longitudinal 
flow of the bulk can change signiflcantly the observed photon energy spectrum via the Doppler 
Effect. Experimental thermal photon data are described by local frame temperature parameter 
which depends on the assumption of longitudinal and radial flow. From the observed hadron elliptic 
flow we estimate the modulation of radial flow parameter versus the angle relative to the reaction 
plane. Based on this we calculate elliptic flow parameter for thermal photons, which was found to 
be very close to that for hadrons. Considering very similar amplitude for hadron elliptic flow at 
LHC and RHIC, we demonstrate that thermal photons at LHC should show large elliptic flow as 
well. All of these considerations are also valid for the low invariant mass dilepton pairs. 

PACS numbers: 25.75.Nq 



Measurements of the photons radiated in nucleus- 
nucleus collisions are of great current interest. This inter- 
est stems from good photon penetration through strongly 
interacting matter. From a theoretical point of view, the 
way to study thermal photons would be to perform a full 
hydrodynamic calculation. Nevertheless, such calcula- 
tions use many assumptions and parameters such as the 
initial energy and density profile, equation of state, treat- 
ment of exact photon production channels, initial tem- 
perature profile, equilibration time, consideration of the 
quark-gluon and hadron phases and many others; these 
assumptions can serve to mask many of the general fea- 
tures of photon emission. One of the common effect on 
top of all these assumptions is the relativistic radial and 
longitudinal expansion of the system. Such expansions 
can induce a Doppler Effect with red and blue shifts for 
radiated thermal photons. Here, we investigate the influ- 
ence of such effects on the emission spectrum for thermal 
photons. We do not pretend to calculate the exact pho- 
ton yield nor specify the actual mechanisms for photon 
production. 

One can distinguish two types of photons: prompt or 
direct photons, produced by primary partons in the ini- 
tial stage of a collision, and thermal photons radiated 
by hot and dense matter after the system thermalizes. 
We limit ourselves to consideration of the latter type. 
We start from a simple system with fixed temperature, 
that expands in the transverse and longitudinal direc- 
tions. Then, we examine how additional assumptions 
can change this simple picture. 

The non-relativistic Doppler Effect for a frequency 
change from a moving source, say, towards an observer is 
w=a;o/(l — v/c), where ujq is the original frequency in 
the rest frame of the source, v is source velocity and c is 
the speed of the signal in the medium. For the relativistic 



case, this formula changes to: 



1 - Pcos0 ■ 



(1) 



Here /3 is speed of the source relative to the speed of light 
and cos6 is the angle between the direction of the moving 
source and the observer. One can see an additional factor 
■\/l — /?2 for the relativistic case. 

For non-relativistic case, if the radiation spectrum has 
a simple Boltzman exponential shape with temperature 
T, 



dN/duJo = exp{—ujQ/T), 



(2) 



using relations loq = uj(l — j3) and dujo/duj = 1-/3, we 
get radiation spectrum seen in the laboratory frame: 

For relativistic kinematics and photons with energy E = 
huj the last formula changes to: 



dN/dE = 



I - PcosO 



Eil- l3cos9) 



)• (4) 



To separate longitudinal and transverse expansion in 
nucleus-nucleus collisions, we describe a boos in the 
transverse direction by velocity /3t or rapidity y^, and 
expansion in longitudinal direction by velocity /3l and 
rapidity y^. The rapidity and velocity are related by the 
formula (3 = tanh{y), where tanh is a hyperbolic tangent. 
The total rapidity y can be calculated via Pifagora's re- 
lation in hyperbolic rapidity space: cosh{y) = cosh{yT) ■ 
cosh{yL). For an observer at mid-rapidity, radial expan- 
sion of the created matter will lead to a blue shift of the 
spectrum when matter moves exactly towards the ob- 
server. In contrast it will induce a red shift for photons 
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produced from the medium, which expands in opposite 
to the observer direction. This expansion in a direction 
opposite to the observer is very similar to astrophysi- 
cal expansion of our Universe. For different cells of the 
medium the vectors for radial flow are oriented at some 
angle to the observer, this effect will be controlled by the 
{StcosO component of the flow. Thus we can write: 



dN/dE 



exp[ 



PtcosO) , 



(5) 



To make quantitive estimates for nucleus-nucleus col- 
lisions we run a simple simulation. We use a density pro- 
files of wounded nucleons for Au-|-Au collisions in x —_u 
plane transverse to the beam direction from work 
which were also used later in paper These profiles 
were generated for different centralities assuming Woods- 
Saxon density distribution of the colliding nuclei. We use 
a linear transverse velocity profile for radial fiow which 
rises from zero at the center of the collision zone to a max- 
imum value of /3r=0.75c This value describes well the 
momentum spectra of identified secondary hadrons pro- 
duced in Au-I-Au collisions at RHIC. The main results 
of our calculation do not change much by changing Pt 
in the range 0.7-0. 8c. Then, each point m x — y plane 
was taken with the weight of relevant participant nucleon 
density. For simplicity we consider the most central 0- 
5% events and set the observer in x direction in reaction 
plane of the collision. For each point we estimate the dis- 
tance from the collision center, calculate the radial boost 
(3{r) — (3t ■ r/R and its projection to the observer direc- 
tion Pt ■ cos9, where the angle 9 in our set is the angle 
relative to the reaction plane (usually marked as (f)). To 
compare with experimental data we have to use the in- 
variant yield by multiplying Eq. [5]by factor E/p^ = 1/E 
for massless photons. 

In Fig. [1] we present results of the calculation. The 
curves were normalized to the first experimental point 
at 1.2 GeV 0. We fit data for ^^=0.5. In this case 
the shape of the experimental spectrum could well be 
described at temperature T=0.17 GeV. There is a sig- 
nificant blue shift of the spectrum at high energy and 
observable red shift at low energy. From the formulas 
presented above, one can see that there is a strong corre- 
lation between parameters T and Pl- Thus, for /3r=0.75c 
the same result could be obtained for different parame- 
ter combinations, like T=0.2 GeV, /3l=0.7c and r=0.3 
GeV, Pl—0-88c or with no longitudinal flow when /3l=0c, 
we get T=0.15 GeV. The experimental data and our es- 
timation for E below 2.5 GeV are very close. Devia- 
tions from experimental points at higher energy could be 
explained by a significant contribution from direct pho- 
tons, whose yield could be estimated by using p-l-p data 
from [3| and scaled by the number of binary nucleon- 
nucleon collisions, thin solid line in figure. The yield 
of direct photons becomes dominant for photon energy 



above 2.7 GeV. In addition, we check how our results 
may change by the varying transverse velocity profile 
(from liner to quadratic) and by introducing a temper- 
ature profile (more hot in the center and colder at the 
edges). We did not find a significant change of the re- 
sults with these assumptions. 

In nucleus-nucleus collisions the produced matter def- 
initely expands in the longitudinal direction and /3l > 0, 
so T=0.15 GeV is the lowest estimate of the mean pos- 
sible temperature of the system. At the same time, it is 
hard to get the actual temperature from the experimental 
spectrum without knowing the value of the longitudinal 
expansion. 
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FIG. 1: Fig.l. Invariant thermal photon yield for central 0- 
5% Au-l-Au collisions versus the observed photon energy in 
laboratory frame for the original source temperature T=0.17 
GeV, /3t=0.75c and /3i=0.5 - thick solid line. The yield was 
normalized to the first experimental point at 1.2 GeV. Points 
are experimental data for central 0-20% Au-|-Au collisions 
from ^ . Thin solid line is for the experimental data for p-|-p 
collisions from the same measurement scaled by the mean 
number of binary nucleon-nucleon collisions at this centrality. 
Dashed line is an estimate of photon yield at very early times 
of the collision just after thermalization occure with temper- 
ature 0.27 GeV but with no flow. 

Another test of our model assumptions is the consider- 
ation of elliptic flow for thermal photons. We can make 
an empirical estimation of the elliptic flow parameter for 
thermal photons using experimental data for hadron el- 
liptic flow. For hadrons produced in A+A collisions we 
can attribute the observed azimuthal anisotropy to an 
amplitude modulation of the radial flow versus the re- 
action plane. Elliptic flow is usually described by the 
parameter V2 as: a{4)) = const • (1 -|- 2 • V2Cos{2(j))), where 
(j> is the angle with respect to the reaction plane. For 
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hadrons,the azimuthal parameter V2-, scales with number 
of constituent quarks and transverse kinetic energy per 
quark § . The value of V2 per quark saturates at about 
W2=0.07 at kinetic energy around 1 GeV per quark for 
mid-central collisions. We use this to estimate the differ- 
ence of radial flow in and out of the reaction plane, /3™ 
and to get a value of V2- The parameter V2 by it- 

self defines the relative variation of the production cross 
section a with respect to the reaction plane orientation. 
So, the difference from the average yield for in-plane or 
out-plane is Aa/a = 2 ■ V2- 

Knowing the shape of the cross section spectrum ver- 
sus the kinetic energy we can estimate what should be 
the difference in kinetic energy to get such cross section 
change. In the kinetic energy range under consideration 
the cross section has an exponential shape with the in- 
verse slope parameter T around 200 MeV for pions and 
300-330 MeV for protons 0, 0. Thus, for the spec- 
trum of constituent quarks the inverse slope parameter 
Tq would be around 100 MeV. If the cross section is de- 
scribed by a simple exponent, a — const ■exp{—Ekin/Tq), 
then its relative variation for an energy change AEkm will 
be Aa/a = AEkm/Tq . To satisfy these changes we get 
the necessary relation with V2 ■ 

AEkrn^2-Tq-V2. (6) 

If a quark with mass nriq moves with velocity /3 and 7 = 
1 / \/l — then its kinetic energy E^in will be E^in = 
niqipi — 1) and we can estimate A7 from the relation 

AEkin = mq ■ Aj = 2 ■ Tq ■ V2. 

At our maximum transverse boost value /3 =0.75 (av- 
eraged over all angles versus the reaction plane), 7=1.51. 
If Tq=0.1 GeV, U2=0.07 and taking a constituent quark 
mass TO = 0.3 GeV we get A7=0.047. This gives es- 
timates of the maximum velocity boost /3™=0.77 and 
/3°"*=0.73. Now we can estimate elliptic flow for thermal 
photons. Fig. [5] shows the value of V2 for 30-35% cen- 
trality Au-|-Au coUisions assuming mq=0.3 GeV (dashed 
line) versus photon energy for the case of /3t=0.75c, 
T=0.2 GeV and /3l=0.7c. We use the similar wounded 
nucleon distribution in the x — y plane as before by cal- 
culating the photon spectrum. For comparison we plot 
the experimental data of V2 for inclusive hadrons ver- 
sus hadron transverse kinetic energy We see signifi- 
cant value of V2 for thermal photons, which is comparable 
with that for hadron V2 and then significantly overshoots 
hadron points at larger energy. We have to take into 
account that in actual experiments, for photon energies 
above 2 GeV, there is a significant contribution from di- 
rect photons (thin line in Fig. [ij for which V2 is about 
zero [7|, this dilutes the total photon V2- We get an es- 
timation of the relative direct photon yield from experi- 
mental results Thus, the measured value of photon V2 
could be surprisingly close (or even larger) to the hadron 
V2- At low energy, where the spectrum is red shifted, V2 



gets negative. Selection of to, = 0.3 GeV is quite arbi- 
trary, so we also made calculations for nig — 0.07 GeV. 
In this case we get A7=0.2, ;3™=0.81 and /3°"*=0.65. 
As one can see, V2 for thermal photons does not change 
much on the assumption of constituent quark mass in 
hadron. Very much the same results we get for combina- 
tions T=0.17 GeV, /3l=0.5c, or T=0.15 GeV and I3l=0. 
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FIG. 2: Fig.2. Azimuthal asymmetry parameter V2 for 
thermal photons versus photon energy in 30-35% centrality 
Au-|-Au collisions. Solid line - for evaluation using the con- 
stituent quark mass = 0.1 GeV, dashed line - for niq = 
0.3 GeV. Experimental points for non-identified hadrons in 
Au+Au collisions at centrality 30-40% versus hadron trans- 
verse kinetic energy are from 5]. Shaded area shows change 
of the result if we take into account the contribution from 
direct photons. 

Resent experimental data for Pb+Pb collisions at LHC 
energy show large hadron V2 with similar values as 
at RHIC. Thus, we can also expect large and at about 
the same value of V2 for thermal photons at LHC and 
red/blue shift effects as well. Similar features, seen for 
thermal photon spectrum, should also be observed for 
low invariant mass dilepton pairs. 

In conclusion, based on a simple assumption that in 
nucleus-nucleus collisions most of the thermal photons 
are produced at the stage when the bulk fiow is already 
developed, we have examined the resulting change in the 
photon spectrum due to longituidinal and transverse ex- 
pansion. The observed picture is very similar to the 
astrophysical expansion. Strong correlation was found 
between the temperature and the longitudinal flow pa- 
rameters. In the absence of longitudinal flow, which is 
definitely not the case, we get the lowest estimate for 
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the thermal photon temperature of about 0.15 GeV. We 
attribute the observed ehiptic flow for hadrons to some 
amplitude modulation of the radial flow. From this we 
estimate the photon flow parameter V2 for mid-central 
collisions, which are found to be large. Because of some 
contribution to the yield from thermal photons, which 
were produced before bulk flow developed, we can con- 
sider our estimates for v2 as upper limits. On the other 
hand, if in experiment we will see significant photon ellip- 
tic flow, it means that the most of photons are produced 
when radial and longitudinal flow are already developed. 
Unfortunately we can not distinguish contributions from 
the quark-gluon or hadron gas stages. 
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